Measurement of a time-periodic magnetic field by rotating coil 
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A novel technique of the measurement data processing is developed which allows to apply the rotating coil method for 
measurement of a dynamic magnetic field, periodic in time. The developed technique allows to obtain time-dependent 
multipoles in a single measurement which takes place during one power cycle of the magnet or the coil rotation period, 
whichever is longer. 
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1. Introduction 

The two-dimensional field harmonic expansion and ro- 
tating coil magnetometer are known for many years and 
are described in many books, reviews and accelerator schools. 
We refer to Jain PQ,[2J and quote a few results which are 
needed for description of the proposed technique of the 
measurement data processing. 

1.1. Plane circular multipoles 

For most practical purposes the magnet field analy- 
sis for accelerator applications can be reduced to a two- 
dimensional problem and the field can be represented by 
the power series expansion (see pQ, sections 1-2): 



B(z) = By + iB x = ^2 c n 
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where B x and B y are the x- and y-components of the mag- 
netic flux density, z — x + iy and reference radius. 
The complex expansion coefficients c„ = b n + ia n in com- 
monly accepted jargon are called "harmonics" or "multi- 
poles" . 

1.2. Measurement of a static field 

The flux $(0) through an arbitrary coil array has the 
form (see [2J, eqn.30) 
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In eqn. ([2| s n are complex "sensitivity factors" , <fi is 
the angle of rotation, "(c.c.)" means complex conjugate. 
We introduced variables 7„ = s n c n which we refer later as 
"reduced multipoles" in order to make further calculations 
less cumbersome. 

In the measurement, flux as a function of the rotation 
angle is evaluated through the time integral of induced 
voltage V(t) (see [2J, eqns. 7,11), readings of the voltage 
integrator are triggered by an angular encoder. In terms 
of time as a function of the rotation angle t(cb) 

tW t(0) t(0) 

/ Vm = J ltt dt = J M = m m (3) 
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As the $(</>) is evaluated in the measurement, "f m and con- 
sequently c m can be found by the Fourier transform: 
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2. Measurement of a time-periodic field 

Let us consider a magnet which is supplied by a peri- 
odic current. If the field is measured by the rotating coil 
method, two cases may be separated, which are considered 
in the following sections. 
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2.1. "Fast coil" case 

We assume the measuring coil rotates uniformly mak- 
ing an integer number of turns M during the magnet cy- 
cling period T. Let us use ^- as the time unit, conse- 
quently the dependence of the "reduced multipole" 7„ on 
time can be represented by a Fourier series: 
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and the flux through the coil vs. time is: 
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The Fourier transform of ^ yields 
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Let us consider a partial sum of Q with n < N, k < K 
and map the a n k to a vector z of length NK as follows: 



z k+{n-l)K — °~nk 
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In the above equation floors'/if) = (j — (j mod K))/K is 
the integer division. The equation fc + nM = m in terms 
of the linear index j defined in (JsJ) can be derived as 



(j mod K)(K - M) + jM = (to - M)K 



(9) 



The smallest solution of the eqn. ([9| is found for m = M, 
j = 0. Let us denote to = m — M. Eqn. ([7| can then 
be rewritten in the form Az = b where the matrix A and 
vector b are defined as: 



(10) 



The case K — M is a trivial one, because A turns to the 
identity matrix, yielding z m = b m . 

In a case K > M, which is of more practical impor- 
tance, the structure of matrix A is more complicated. Let 
us consider as a case study A for K = 3, M = 1 and 
N = 3: 
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One can see the matrix contains K N — (K — M) (N — 1) = 
K + M(N — 1) rows. It is filled with a regular pattern 



made of N parallel to the diagonal groups of ones, each 
group of length K. Positions of ones in n th group are 
given by equation to = j — (n — 1)M. Each column in the 
matrix contains exactly one "1" , which means that each 
unknown enters only into one equation, thus all equations 
are decoupled. 

As the number of equations (rows) is less than the 
number of unknowns, the system must be supplemented 
by more equations to have a unique solution. These ad- 
ditional equations can be obtained using the translation 
properties of a n k and b m . 

Each equation with L > 1 terms has, after mapping 
back from Zj to a n k, the following form: 
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In equation ( 12 1 j m j n is the minimal index with z- mapping 



index n and the maximal value of index k in the corre- 
sponding a n k set, respectively. 

If the measurement clock is shifted by r, from ^ and 



( 10 ) it follows that o~ n k and b rn are transformed as: 
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In terms of the equation for Zj , from the original equation 
( 12 ) we can construct another one: 
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Let us construct L—l additional equations as in ( 14 ) with 
different r; with / < L. Denoting pi = er lMTl and setting 
Po = 1 we see that the matrix V built from rows of coef- 
ficients of each equation is a square Vandermonde matrix 
[3] with a known determinant: 
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0<i<j<L 

If all pi arc chosen to be different (a natural choice: pi = 
e~ 27r4 r) 7 V has a non-zero determinant, therefore it has 
an inverse, which for the Vandermonde matrix can be ex- 
plicitly expressed @]. 

Applying the above technique to the case study we ob- 
tain the following set of equations: 
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As all matrices in ( 16 ) are invertible, the case study is 
solved. This technique can be applied for any N, M and 
K > M. 

Instead of inverting the NK x NK matrix, the problem 
is reduced to inverting of a number of matrices of smaller 
dimensions. It can be shown for the general case, that the 
maximum matrix dimension does not exceed (K — M + 
1) x (K - M + 1). 

2.2. "Fast magnet" case 

Let us assume the measuring coil rotates uniformly and 
the magnet is cycled M times during one turn of the coil. 
In this case we denote as T the coil rotating period and 
measure time in units of . For the purpose of a partial 



reuse of results obtained in section 2.1 let us swap the 
meaning of indices: now k stands for the multipole number 
and n for term number in the Fourier expansion of the 
multipole: 
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The flux through the coil vs. time is then given by equation 
which is similar to Q, except for the low limits of sums: 
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The Fourier transform of ( 18 1 yields 

$ e~ mt dt = ^2^2 S rh : (k+Mn)Vkn (19) 
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Let us consider a partial sum of (18) with n < N , k < K 



and map the <j n k to a vector z of length NK as follows: 

Zk-l+nK = Okn ^ Zj = & j mod K+1< floor(-i) , . 

(20) 

< j < NK 

The equation k + nM — rh in terms of the linear index j 



defined in ( 20 ) can be derived as 



( j mod K) (K - M) + jM ={m-l)K (21) 



The smallest solution of the cqn. (21 1 is found for m = 1, 
j = 0. Let us denote m = rh— 1, then ( 19 1 can be rewritten 
in the form Az = b where the matrix A and vector b are 
defined as: 
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Comparing ( 22 1 and ( 10 ) one can see that the same matrix 
A is obtained (albeit after swapping meanings of k and 
n). Therefore the consideration of the "fast magnet" case 
is very similar to the one of "fast coil" . 

The case K = M is a trivial one, because A turns to 
the identity matrix, yielding z m = b m . 

For K > M each equation with L > 1 terms has, after 
mapping from Zj to (Jkm the following form: 
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If the measurement clock is shifted by t, a^ n and b m are 
transformed as: 
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In terms of equation for Zj , from the original equation ( 23 1 
we can construct another one: 
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Consequently, exactly as for the "fast coil" case, from each 
equation with L > 1 terms we construct L — 1 additional 
equations using the Ukn and b m time translation proper- 
ties, obtaining a linear system with a Vandermonde ma- 
trix, which has a non-zero determinant. 

Just as in the "fast coil" case, the problem is reduced 
to inverting a number of matrices with dimensions not 
exceeding (K - M + 1) x (K - M + 1). 

3. Discussion 



For the "fast magnet" case (Section 2.2) a good syn- 
chronization of the coil rotation and the magnet cycling is 
essential. As a possible implementation one can consider 
the magnet power supply driven by a waveform generator 
which is clocked by the angular encoder pulses. 

Much looser requirements can be applied for the "fast 
The magnet cycle starts from 
If multipoles 



coil" case (Section 2.1). 
the injection plateau and comes back to it. 
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stay constant at the injection plateau then expulsion of 
a part of the plateau, followed by "gluing" of remaining 
parts preserves continuity of each multipoles as a function 
of time. 

Therefore in the "fast coil" case the measurement can 
be started at any point shortly before the ramp up, and 
ended at any point during the injection plateau after ramp 
down when the coil comes back to the start angle. In other 
words, the magnet cycle is not required to be a multiple 
of the coil rotation period. This option, to choose a "ficti- 
tious" magnet cycle, may be exploited. It may be shown 
for a cycle with the same ramp up and ramp down rates, 
that the closer the "fictitious" injection plateau duration 
is to the duration of the cycle flat-top the lesser are higher 
harmonics in the Fourier expansion 

The developed data processing technique operates with 
truncated expansions (partial sums) to represent the field 
and multipoles. To check if upper limits of sums are large 
enough to provide the precision required they may be var- 
ied and the effect of variations on obtained results must 
be analyzed. 



4. Conclusions 

The data processing technique has been developed which 
allows the use of the rotating coil method for the measure- 
ment of a dynamic time-periodic field. This technique al- 
lows to obtain time-dependent multipoles in a single mea- 
surement which lasts one power cycle of the magnet or the 
coil rotation period, whichever is longer. 
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